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Generalized hydrodynamic theory, which does not rest on the requirement of a local equilibrium, 
is derived in the long-wave limit of a kinetic equation. The theory bridges the whole frequency range 
between the quasistatic (Navier-Stokes) hydrodynamics and the high frequency (Vlasov) collisionless 
limit. In addition to pressure and velocity the theory includes new macroscopic tensor variables. In a 
linear approximation these variables describe an effective shear stress of a liquid and the generalized 
hydrodynamics recovers the Maxwellian theory of highly viscous fluids - the media behaving as solids 
on a short time scale, but as viscous fluids on long time intervals. It is shown that the generalized 
hydrodynamics can be applied to the Landau theory of Fermi liquid. Illustrative results for collective 
modes in confined systems are given, which show that nonequilibrium effects qualitatively change 
the collective dynamics in comparison with the predictions of the heuristic Bloch's hydrodynamics. 
Earlier improvements of the Bloch theory are critically reconsidered. 
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I. INTRODUCTION 

The idea to apply a macroscopic hydrodynamical de- 
scription to the collective dynamics of the inhomogeneous 
many-electron systems has been first suggested in a fa- 
mous work by Bloch in 1933.1 In this paper he intro- 
duced a hydrodynamic theory of a degenerate Fermi gas 
as a simplest phenomenological extension of the Thomas- 
Fermi model to dynamical regime. Since application of 
nonequilibrium many-body theory to inhomogeneous sys- 
tems is extremely complex, the Bloch's hydrodynamic 
theory (BHT) remains popular and is in use until now. 
Since 1933 theJBiHT has been applied to variety of ki- 
netic problems pu albeit due to the heuristic nature of 
the BHT, its relation to the kinetic theory and the range 
of applicability remained unclear. 

A common sense of hydrodynamical theory presumes 
an existence of a local statistical equilibrium in every 
point of space - the condition which allows to express the 
kinetic properties of a many-body system in simple-terms 
of pressure, density and temperature distributions. U This 
condition is, however, rarely met in electron gas where 
the electron-electron collisions commonly play a minor 
role. This circumstance is ignored in the BHT. Al- 
though, the-, applications of the Bloch's hydrodynamical 
theoryBiaiJ have shown its usefulness, its failures are 
also well known, irupactiewilar a prediction of a wrong 
plasmon dispersion .q~TD , E£H13 

Is it possible to derive a macroscopic hydrodynamics- 
like theory which would be valid beyond the condition 
of a local equilibrium? What is the role of collisions and 
how is the static (i.e. collision dominated) limit as well as 
the high frequency (collisionless) limit recovered? These 
are the questionS|-s«e address in the present paper. In the 
short publicationEj we showed that such hydrodynamics 
existed. In the paper Ref. [1^, for the sake of simplicity 
we carried out the derivation of the hydrodynamics equa- 



tions only in one dimensional case. We also restricted the 
theory to the second order in the long-wave limit of a ki- 
netic equation, which left out the thermal conductivity 
and consequently the thermal equilibration processes. 

In the present paper we give a new general derivation 
and extend the theory to include the heat transport. We 
show that the generalized hydrodynamics can be based 
on Landau theory of Fermi liquid and relate macroscopic 
parameters in hydrodynamics equations (i.e. shear and 
bulk moduli), to microscopic Landau parameters. We 
also present transparent analitical examples of the calcu- 
lation of collective modes in confined systems of different 
dimensionality. 

Let us first consider the mathematical structure of 
BHT and the formal microscopic restrictions of its ap- 
plicability. The Bloch's original idea was to use the set 
of hydrodynamics equations of an ideal charged liquid to 
describe the dynamics of a Fermi gas. Only the macro- 
scopic variables - electron density n(r, t), velocity v(r, t), 
pressure P, and electrostatic potential (p(r, t), enter the 
equations: 

D t n + nVv = 0, (1) 
mnD t v + VP — enVtp = 0, (2) 
VV = 47ren - Airp ext , (3) 

where p ex t is an external charge and D t — dt + vV is a 
substantial derivative. This set of equations (continuity 
equation ([!]), Euler (||) and Poisson (||) equations) be- 
comes complete when the equation of state is added. In 
the original paperlll Bloch identified P with the kinetic 
pressure of a degenerate Fermi gas: 

P=— (37r 2 ) 2 / 3 n 5 / 3 (4) 
5m 

Further improvements of BHT addressed only the equa- 
tion of state (inclusion of eschange, correlation and quan- 
tum gradient corrections). □ 
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Sometimes BHT is viewed as an approximate extension, 
of the density functional theory to the dynamic regime.u 
However, as we have mentioned above, the hydrodynam- 
ics theory (|l])-(^|) exhibits an inconsistency which-Jiaa 
been reflected in many papers and textbooks fll3'll2rtHl 
The system of equations (0)-(§) gives the wrong veloc- 
ity coefficient v 2 , in the plasmon dispersion law lu 2 — 
u)p + v 2 q 2 . Instead of the correct result |i;f, (vp is the 
Fermi velocity) of the linear response theory (RPA) the 
BHT leads to the value v 2 = \v"p. This inconsistency 
does not depend on a degree of degeneracy of an electron 
gas and is the general property of the hydrodynamics the- 
ory which is based on the equilibrium equation of state. 
For example, at arbitrary degeneracy the result of hydro- 
dynamics is Vq = Vg, where v s is the velocity of sound, 
whereas in the linear response theory Vq eauals to the 
mean square of the particle velocity < v 2 >.E3 

It has been realizedE3 that this discrepancy originates 
from the assumption of a local equilibrium, which under- 
lies the common hydrodynamic theory!] The assumption 
allows to reduce the kinetic equation for the distribution 
function / p (r, t) 

dtf P + -V/ p + = J p [/ p ] (5) 

m op 

to a set of equations for macroscopic variables 
n(r, t),v(r, t) and, in general case, temperature T(r,t). 
The requirement of the local equilibrium is fulfilled if 
the characteristic time of the process t pr ~ l/u> is much 
longer than the inverse collision frequency l/v c and the 
typical length of the inhomogeneity L is greater than 
the meap-ifree path I ~ u/v c (u is the average particle 
velocity 

uu/v c < 1, u/Lv c < 1. (6) 

In a zero order with respect to parameters (^) Eq. (|^) 
reduces to I p [f p ] — 0, which means that the distribution 
function / p and, consequently, the equation of state has 
a locally equilibrium form. As a result, the equations for 
the first three moments of the distribution function i.e. 
density n, current j = nv and stress tensor i^y = SijP 
(which is diagonal in the local equilibrium) constitute the 
closed set of hydrodynamics equations for an ideal liquid. 

It is clear that from the microscopic point of view BHT 
cannot be a consistent theory since it extends the colli- 
sion dominated hydrodynamics (y c — > oo) to the electron 
gas where the collisionless (Vlasov) limit (y c — > 0) is most 
common. Due to the high frequency of the plasma waves 
at least the first of inequalities (g) is strongly violated 
and the tensor structure of as well as that of higher 
moments becomes important. In a degenerate case these 
tensors describe deformation of the Fermi sphere - the 
effect completely ignored in the BHT. It has been shown 
recentlylij that within a linear response theory this effect 
leads to an effective shear modulus of a liquid. 

In this paper we present a generalized hydrodynam- 
ics which remains valid far beyond the local equilibrium, 



when both conditions (|6|) are strongly violated. The the- 
ory is restricted to the long- wave limit of the kinetic equa- 
tion and requires inclusion of the tensors of higher mo- 
ments. In Sec. II we derive a general hierarchy of equa- 
tions for the moments of the distribution function and 
introduce a regular procedure of truncation of the infinite 
chain of equations, which does not require the assump- 
tion of a local equilibrium. We also discuss the physical 
meaning of restrictions which make possible the trunca- 
tion. In Sec. Ill we construct a hydrodynamic theory of 
the second order in the long wavelength expansion and 
show a relation of the generalized hydrodynamics to the 
theory of elasticity and to the theory of highly viscous 
fluids. Sec. IV is devoted to the generalized hydrody- 
namics of a charged Fermi-liquid on the basis of Landau 
theory. We find the relationship between macroscopic hy- 
drodynamical parameters and microscopic Fermi-liquid 
parameters. Sec. V contains examples of application 
of the generalized hydrodynamics to collective oscilla- 
tions of confined systems. The purpose of this section 
is purely illustrative. To show clearly the nonequilib- 
rium effects we select an analytically solvable model of a 
charged Fermi liquid in a harmonic potential of different 
dimensionality. We show that deviations from the local 
equilibrium strongly shift the excitation frequencies and 
substantially change the structure of collective modes. In 
Sec. VI we extend the theory to the fourth order in the 
long wavelength expansion. It is shown that in the low 
frequency limit the fourth-order theory transforms into 
the set of equations of a classical hydrodynamics, which 
includes viscosity and a heat transport contribution. In 
Sec. VII we summarize our results. 



II. THE HIERARCHY OF EQUATIONS FOR THE 
MOMENTS 

The kinetic equation (^|) can be transformed in an in- 
finite chain of-e/quations for the moments of the distribu- 
tion function. til The zeroth and the first moment are the 
particle density n(r,t) and the current density j(r, t): 

n(r, t)=J2 /p(r, *), *) = £ £ / P (r, *) (?) 

p p 

The velocity field can be defined in a usual way as 
v = j/rt. In general, the moment of the fcth order is 
a tensor of fcth rank 

p 

It is a common practice to separate the macroscopic mo- 
tion of the liquid from the relative motion of particles via 
transformation to the Lagrange comoving frame. In this 
frame the macroscopic velocity is equal to zero and the 
particle density n(r, t) is the same as in the laboratory 
frame. The transformed distribution function f£ (r, t) 
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is related to the distribution function in the laboratory 
frame / p (r,i) as follows 

fp = fp+mv 

and apparently has the properties 

= £/ p (r,t) = n(r,t), (9) 
p p 

E £ / P L ( r < *) = E £ /p( r ' *) - v ( r ' *M r > *) = °- ( 10 ) 

Z ^ 771 * Z ^ 777 



Eqs. (|9j) and (|10|) can be considered as a definition of the 
comoving frame. 

The higher moments of the transformed distribution 
function are defined similar to (|J): 



(fc) _ r(fc) 



(fc) 

Equations for moments i can be easily derived from 
the kinetic equation for the distribution function f£: 

df L 

-(mDtv-eVrt-jg (12) 

where D t = dt + vV is a substantial derivative. The 
zcroth and the first moments of this equation give re- 
spectively the continuity equation and the equation for 
the velocity v: 



D t n + nViVi — 

(2) 

mnD t Vi + VjL- — enVif = 



(13) 
(14) 



— ^p ii ...p ikPj —V j v i = 



In the derivation of Eqs. (|l3|)-(|lj) the constraints due 
to conservation of particles, momentum and energy in a 
collision process 



E^ P = o, e^ 7 p = ' Ep 2/ p = ' 



(15) 



(2) 

have been used. The second moment L;„- is a stress ten- 
sor, which is related to the pressure P of an electron gas 

(2) 

as TrL^j = 3P. This tensor can be decomposed into a 
scalar and a traceless parts 



L 



(2) 



PS* 



Tr-Ki 



0. 



(16) 



The equation of motion (hj) allows to eliminate the elec- 



trostatic potential tp from the kinetic equation (12) 



Dtf, 



m 



-Vifp 



PjVjVi 



V.,L 



(2)' 



3 i] 



dpi 



= Ip (17) 



and, consequently, from all equations for higher mo- 
ments. The term in brackets in Eq. ([Tt]) has a clear 
physical meaning. This is a force which acts on a parti- 
cle in the Lagrange frame. 

Calculation of the fcth moment of the kinetic equation 
(O) leads to the equation for the fcth moment of the dis- 



tribution function L 



a-) 



The fcth moment of the first 



i — □ 

term in the left hand side of Eq. ( |17| ) gives the substantial 
derivative of the fcth moment of the distribution function 
L/ fe ). Similarly the divergency of the fc + 1th moment 
L(fc+i) comes f r om the second term. Contributions from 
the last two terms on the left have more complicated ten- 
sor structure. For example, the fcth moment of the first 
term in the brackets in Eq. (HJ) is 



dpi 



—y 

,fc-i 



dpi 



^in L i 2 ...i k j 



^2 L i 1 i 3 ...i k j 



'V'< 



(18) 



The contribution from the last term in the left-hand 
side of Eq. (|l7j) has a similar structure but contains the 
fc- 1th moment L^"" 1 ) instead of L/ fc ) in Eq. (18). Since 



the first moment of the distribution function f£ equals to 
zero (see Eq.(fl0|)) this contribution vanishes in equation 
for the second moment. Hence, in the case of fc = 2 one 
has the following equation for the moments: 



(2) 

where 1^ is the second moment of the collision integral. 
By the definition (|ll|), all moments L\ i are sym- 



(2) 



Lfv k v k +L%V kVj 



'(2) 



,(3) 
J ijk 



r(2) 



(19) 



I 

metric with respect to transmutation of all indexes. All 
equations for moments must preserve this symmetry. To 
explicitly incorporate this condition it is convenient to 
introduce special notations for two possible symmetric 

(k) 

products of a symmetric fc-th rank tensor • and a 
vector cij. The first one is the dot product, which de- 
creases the rank of a tensor from fc to fc — 1 and corre- 
sponds to the contraction of the vector index j with one 
of tensor indexes 



T (fc) 



T, 



(fc) 



(20) 



We also define a symmetric direct product of and 
a which gives a tensor of the rank fc + 1 and is equal to 
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the sum over all transmutations of the vector index j and 
tensor indexes i\...ik 



Lmax{ui, u c } 



< 1. 



(26) 



T {k) a- 



T iiL k _ d "i k - (21) 



These notations allow us to rewrite the moments of dif- 
ferent terms of the kinetic equation (|l7| ) in a compact 
way. For example, the A;th moment of the first term in 
the round brackets in Eq. (|l7|), which is given by Eq. 
(fLq) above, can be expressed as 



— ^p ii ... PikPj —V j v i 



dpi 

= lW(V-v)+{(lW-V)®v} s . 

Similarly, the moment of the second term in the brackets 
in Eq. @ is 

^E«.-«.fv,Lg' = i{L<-.«(V.L«) }s 



Finally the system of equations for the moments takes 
the form 



D t n + nV-v = 0, 



(22) 
(23) 



mnD t v + V ■ L (2) - en\7ip = 0, 
AL (2) + L( 2 )(V • v) + {(l/ 2 > ■ V) ® v} s 

+ V • L< 3 > = I (2) (24) 
AL (fc) + L w (V • v) + {(L (fc) • V) ® v} s + V • L (fe+1) 

_^{L( fc - 1 )®(V-L( 2 ))} s =lW, (25) 



mn 



where the last equation ( pq ) is valid for k > 2. 

In a collision dominated case when conditions (^) are 
fulfilled, the high moments are small and the infinite 
chain of equations (p2[)-(p5|) can be truncated (Chapman- 
Enskog or Grad methods, see Ref. As the traceless 
part 7Tjj of the stress tensor (Q) and the third moment 



(3) 

tensor L\- k vanish for the locally equilibrium distribution 
function, they remain small under conditions (H) . In this 

(3) 

case tensors %ij and L\-' k describe, respectively, viscosity 
and thermal conductivity in a collision dominated liquid. 
However, if conditions (||) are violated, itij and L^ k as 
well as the higher-order moments are not small. Yet the 
infinite chain of equations can be decoupled if all physical 
quantities are slow varying functions of r. It is seen from 
Eq. ( p^ ) that the third moment L^X enters only under a 
spatial derivative, hence its contribution is proportional 
to 1/L. This remains true for all higher moments in 
higher-order equations (^5|) . Thus the truncation of the 
chain can be guaranteed by the smallness of the gradients 
of the moments instead of the smallness of the moments 
itself. The dimensionless parameter of this expansion is 



It is clear that both inequalities (||) can be violated while 
the condition (^6|) is fulfilled. For collisionless nonde- 
generate plasma the possibility to decouple the chain of 
equations for moments based on the condition similar to 
Eq. (p6| ) was first pointed out in Ref. 16, where it was 
referred to as "the low temperature approximation" . 

To justify the consistency of this truncation procedure 
we estimate the order of magnitude of the terms in Eq. 
(^||) . According to Eq. (^2[) , the spatial derivative of the 
velocity field v is of the order of the inverse characteristic 
time l/tp r ~ u> hence the first three terms in Eq. fl25| ) 
are proportional to to. The right-hand side of this equa- 
tion is evidently proportional to the collision frequency 
v c . The last two terms in the left-hand side contain only 
spatial derivatives and are of the order of u/L. Thus 
the contribution of Lw in the equation for L/ fc_1 ) con- 
tains an additional smallness 7 (p6t), which means that 
the correction originating from L™ to the equation for 
the density (|2^) is of the order "f k . Therefore to obtain 
a theory valid up to 7 fc one should keep k + 1 equations 
for L/ 9 ' (0 < q < k) with the contribution of the k + 1-th 
moment and the spatial derivative of I/ 2 ' being omitted 
in the last equation. The resulting theory is a gener- 
alization of hydrodynamics which is valid far from the 
equilibrium. 

It should be mentioned that, in principle, the described 
truncation procedure may not provide the closed system 
of equations because the moments of the collision integral 
are, in general, the functionals of the distribution func- 
tion. However, for a high frequency process (<jj/v c 1) 
the contribution from the collision term disappears and 
we obtain the closed set of equations of the " collisionless 
hydrodynamics" which corresponds to the Vlasov limit of 
the kinetic equation. The collision terms become eventu- 
ally important on the long time scale, when the system 
approaches the equilibrium. In this limit the moments 
of the collision integral are linear functions of the corre- 
sponding moments of the distribution function.lifl Hence 
under the condition (uj/i> c <C 1) they can be presented 
in a close form. Matching these two limits gives the hy- 
drodynamical theory which correctly describes the high- 
and the low-frequency limit and is approximately valid 
for all values of lo/v c . 

There is an obvious relation of the generalized hydro- 
dynamics to the standard linear response theory. In the 
high-frequency region the collisionless hydrodynamics is 
a regular expansion of Vlasov equation in terms of param- 
eter (p(|). Consequently in the case of a weak perturba- 
tion the linearized hydrodynamics equations should give 
the same results as the linearized Vlasov equation. For 
example, the response function obtained from the lin- 
earized system (^)-(|2^) which is truncated at /cth equa- 
tion should coincide with the microscopic response func- 
tion up to {u/Loj) k ~ (uq/ut) k (where q ~ 1/L is a wave 
vector of a perturbation). However the hydrodynamics in 
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its general form goes beyond the linear response theory 
(RPA) since it is valid for nonlinear regime. 

In this section we introduced parameter 7, which gov- 
erns the decoupling procedure. The smallness of this pa- 
rameter restricts the region of applicability of the gen- 
eralized hydrodynamics. The condition ( p6| ) physically 
means that a particle passes the characteristic length 
scale L in a time much longer than the typical time of 
the process ~ 1/uj or a collisional time ~ \/u c . In this 
case there are no particles which move in resonance with 
a collective motion and may provide an energy exchange 
between single-particle and collective excitation. In other 
words, the contribution of the Landau damping to the 
evolution of the system is small. This is the physical 
restriction of the applicability of the generalized hydro- 
dynamics which is (as any hydrodynamics) the theory of 
a collective motion. 



III. THE SECOND-ORDER THEORY. 
HYDRODYNAMICS OF A DEGENERATE 
FERMI GAS 

The theory of the second order with respect to param- 
eter 7 ( |26|) corresponds to neglect the third moment L/ 3 ' 
in Eq. (J24Q. Let us decompose the stress tensor into the 
scalar and the traceless parts (|l^) and write the system 
of equations in terms of the fields P and tt: 



D t n + nV • v = 0, 

to D t v + VP + V ■ 7r - eriVip = 0, 

D t P + |pV.v + |(ir.V).v = 0, 



(27) 
(28) 



(29) 



D t TT + 7r(V • V) + {(tT • V) ® V} S - 1-(7T ■ V) ■ V 



P ( {V®v} s - l^V • v ) = I (2 '> 



(30) 



where 1 is the unit tensor. Eqs. (2£) and ( pOj ) are the 
trace and the traceless part of Eq. ( 24 ) respectively. Due 
to the conservation of energy in a collision process the 
second moment of the collision integral 1^ is a traceless 
tensor. Consequently it contributes only to Eq. (30). 
Eqs. (p7|)-(p0|) are the hydrodynamics equations which 
contain two scalars n and P, one vector v and one trace- 
less tensor tt. For a charged liquid one has to add the 
Poisson equation (||) for the scalar potential ip. 

Eigenmodes of a charge liquid (plasma waves) lie in 
a high frequency region where the collisionless limit is 
reached. In Sec. IV (see also Ref. |l3|) it is shown that 
the system of equation (p7|)-(|3C|) leads to the correct plas- 
mon dispersion up to the second order of parameter 7. 
This is a reflection of the fact that the linearized Eqs. 
(f27j)-(|30|) correspond to a correct linear response theory 
up to the 7 2 . 

In the low frequency limit an expression for collision 
integral is needed. To demonstrate how a long-time scale 



behavior is recovered, we take I p in a Krook-Bhatnager- 
Gross (KBG) approximation (see, i.e. Ref. [n]) 



/p[/p] = -M/ P (r,t)-/p(r,i)), 



(31) 



where /p 7 is a local equilibrium Fermi function with 
position-dependent velocity, chemical potential /x(r, t) 
and temperature T(r,i). By the definition the function 
fp is chosen to give the same values of velocity v(r, t), 
density n(r,t) and pressure P(r,t), as the exact distri- 
bution function f p (r,t): 



n(r,t)=X;/p=E/i 
p p 

j(r,t)=x; £ /p = E £ /p 

^— ' TO TO ^ 

p p 

^m)=e£w=e£/ p * 



3to 



F 

3m' p+mv 



(32) 
(33) 
(34) 



Eqs. ( |32[|34| ) guarantee that the collision integral ( |3l| ) 
satisfies the general propertie s (|l5| ). 

The second moment of J p fl31|) is equal to 



1^ = -u c (l^-1P f (u,T)), 



(35) 



where Pf{it,,T) is the pressure of a Fermi gas with the 
distribution function f£. By definition (|34|), this pressure 

equals to the exact pressure Pp(n,T) = P = ^Trl/ 2 ), 
hence the expression in the brackets in ( Pq) e quals to the 
traceless tensor tt. Consequently, Eq. (p5ty for I^ 2 ) is 
equivalent to two equations 



I< 2 > = -^cTT, 

P=P F (n,T). 



(36) 
(37) 



Since Pf(w, T) is an equilibrium pressure for given n and 
T, the equation ( |37| ) may seem to be related to the as- 
sumption of the local equilibrium. In fact, it is not. Eq. 
( |37| ) simply introduces the new independent scalar vari- 
able T(r,t) instead of P(r,t). The system can still be in 
an arbitrary (up to 7 2 ) nonequilibrium state. In the sec- 
ond order theory such state is uniquely described by two 
scalar, one vector and one traceless tensor functions. The 
traceless second rank tensor tt is responsible for devia- 
tions from the local equilibrium. In a degenerate Fermi 
system this tensor describes deviations of the shape of 
the Fermi surface from sphere. 

The system of equations (^7| )-(^0|) , with the collision 
term (|3^) and P from Eq. (|37|) transforms into the com- 
mon hydrodynamics theory in the limit ui/v c <C 1. In- 
deed, as we already mentioned above the spatial deriva- 
tive of v has the order of magnitude of to. Hence, the 
first four terms in the left hand side of (^Tj) are pro- 
portional to W7T, the last term in the left hand side is 
proportional to ujP and the right hand side is equal to 
v c tt. Consequently, in a zero order of lu/v c Eq. (|30|) 
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leads to 7Tjj = 0, and Eqs.(|27D-(|29[) become identical to 
hydrodynamics equations for an ideal liquid: 



D t n + nVv = 0, 
mnD t v + VP — enVip = 0, 

D 4 P+-PV-v = 0. 



(38) 
(39) 

(40) 



Here Eq. (|39|) is the Euler equation and Eq. ([40[) is the 
equation for the conservation of energy. For a degener- 
ate Fermi gas this set of equations exactly corresponds 
to the BHT. We emphasize again, that it is valid only in 
an extremely collision dominated regime. 

In the first order of u>/v c <C 1 one has to neglect the 
first four terms in the left hand side of Eq. (^) in com- 
parison with the last term which is proportional to the 
pressure P. Thus, the first order solution of Eq. ( J30| ) 
takes the form of the viscosity tensor 



P 



(41) 



with the viscosity coefficient ?/ = P/v c . In this case Eqs. 
( p8| ) and (|2|) are equivalent to Navier-Stokes equation 
and the equation for energy conservation in a viscous liq- 
uid respectively. In Eq. ( f4l|) we have turned back to 
common tensor notations to make the relation to the or- 
dinary viscous hydrodynamics more transparent. 

The formulas ( ^6|) , ( |37| ) were obtained using the KBG 
collision integral (pa). They are, however, more general 
than the KBG approximation itself. As it has been men- 
tioned at the end of the previous section, the reason is 
that in the low frequency range, where collisions are im- 

(2) 

portant, the second moment P- is always a linear func- 



tion Of 7T, 



.0 



The coefficient v r can thus be considered as 



a phenomenological parameter, which is related to vis- 
cosity r\ as v c — P/rj. 

Eqs. ©-© with l( 2 ) from @ and P from @ con- 



stitute the closed set of equations of generalized hydro- 
dynamics in the second order of 7. This set of equations 
gives a correct description of the high-frequency collision- 
less regime. It also leads to the classical hydrodynamics 
of a viscous fluid in a low- frequency limit. Consequently 
it should be valid, with a reasonable accuracy, for inter- 
mediate regime. 

However, the second order theory has an inherent in- 
consistency. Though in a low frequency limit the viscous 
term is recovered correctly, the thermal conductivity con- 
tribution V 2 T (Ref. ^|) is still missing in the energy con- 
servation equation. Without this contribution the correct 
static limit VT(r) = cannot be recovered. In fact, the 
term V 2 T corresponds to correction of the fourth order 
of 7 in the continuity equation. Physically this means an 
absence of a dissipative flow of energy with an accuracy 
up to 7 2 . To include the thermal conductivity one has to 
consider the third and the fourth moments and neglect 
the fifth moment (see Sec. VI bellow). 



The theory of the second order in 7 is obviously ap- 
plicable for a description of dynamics of a dense degen- 
erate Fermi gas. In this case the condition T/Ep <C 1 
(Ep is the local Fermi energy) is always satisfied and 
it is unnecessary to take into account the third and 
the fourth moments. At T = Eq. (|37|) reduces to 
P = ^-(3n 2 ) 2 / 3 n 5 / 3 and Eq. @ transforms into 



V) 



0. 



The equation for pressure P and the condition 
gether with the system of differential equations 



D t n + nV • v = 0, 
mD t v + VP + V ■ 7T 



enVip = 0, 
D t n + tt(V • v) + {(tt • V) <g> v} s - 



ifv-v 



l-(TT-V) 



P {V <g> v}< 



(42) 
|) to- 

(43) 
(44) 

v 

(45) 



provide a complete system of equations of the generalized 
hydrodynamics for a degenerate Fermi gas. 

Most applications of the Bloch's hydrodynamics con- 
cerned just the case of the degenerate Fermi gas. The 
set of equations ((42|)-(f4"5|) of generalized hydrodynamics 
should replace the BHT, which is incorrect except the 
collision dominated limit or the case of spatially homo- 
geneous system. 

There is an interesting relation of the linearized version 
of the generalized hydrodynamics in a high-frequency 
(collisionless) limit to the theory of elasticity. In the limit 
u>/v c ^>l linearization of the system (p7^-(|30|) gives 



d t Sn + n VkVk = 
mnodtVi + Vi<5P + V^TTy — engViip = 

d t SP+^PoV k v k = 

/ 2 
dtWij + Pq Vj^j + VjVi - -5 i:j V k v k 



0, 



(46) 
(47) 

(48) 
(49) 



where Sn and SP are deviations from the equilibrium 
density hq and pressure Pq. Let us introduce a dis- 
placement vector u(r, t) as Sn = — rioVu. The continu- 
ity equation ( [46|) gives a usual relation between velocity 
and displacement dtu = v. Introducing a stress tensor 



-SPSi 



nij one can rewrite Eq. (47) as 



mnod^Ui — X7j<Ji 



en Vitp = 0. 



(50) 



The relationship of the tensor cry to the displacement itj 
follows from Eqs. ([pL |49|) and takes the same form as in 
the elasticity theor} 



er i7 - = KV k u k 5i 



^jVfcMfe ) , (51) 



where the bulk modulus K and the shear modulus (i of 
an electron gas are K = |Pq and \i = P . Physically the 



G 



bulk modulus is responsible for the increase of the en- 
ergy which is caused by the local change of the occupied 
volume in a momentum space, whereas the shear stress 
describes the deviation of the shape of this volume from a 
sphere. In the next section we derive the hydrodynamics 
for a Fermi liquid and find the correlation contribution to 
bulk and shear modulus. The elastic description of Fermi 
systems has been recently discussed in Ref. [l5| within the 
linear response theory. 

Inclusion of collisions violates the exact correspon- 
dence to the elasticity theory. In this case the shear stress 
tensor cry — ^5ijO~kk = ^^ij should be determined from 
the equation 



d t iTij + (J-dt ( Vjitj + VjM, - -Si 3 \7 k Uk 



(52) 

which follows from Eqs. ( fl9| ) and ([36|). Eq. ( p2( ) exactly 
coincides with the equation for the shear stress tensor in 
the phenpsnenological theory of highly viscous fluids by 
MaxwellJIa These fluids (for instance, glycerin or resin) 
behave as solids at short intervals of time, but as viscous 
liquids on a large time scale (see Ref. |l8|). The general- 
ized hydrodynamics thus provide a surprising but clear 
analogy between an electron gas and Maxwellian highly 
viscous fluids. 



IV. THE GENERALIZED HYDRODYNAMICS OF 
A FERMI LIQUID 



to the comoving frame. The kinetic equation for the dis- 
tribution function in this frame looks similar to Eq. 
(|l2|), but with an additional correlation force: 



T d£ » T , s 9n p 



(mAv - eW(p + V£p) 



dnh, T 

P_ jL 

dp ~ p 



(54) 



According to the Landau assumption, the quasiparticle 
energy is a linear functional of the deviation of the dis- 
tribution function from the Fermi step function: 



e L = e<°> 



■£/i 

p' 

-5nl 



n^=6(p- PF ), 



(55) 
(56) 



where £ p ^ is the quasiparticle energy in the undisturbed 
system and f p / p is the quasiparticle-quasiparticlc inter- 
action. Linearization of the kinetic equation (54) leads 
to the Landau-Silin equationjlia written in a comoving 
frame 



d t Sn^ + 



dn 



(0) 



dpi 



- ( md t Vi - eVi^p + fp'pViSn^, 



dn 



(0) 



dpi 



= lt (57) 



The quasiparticle mass to*, which enters equation (|5 
is related to the bare fermion mass to as follows^ 



In this section we develop the hydrodynamic descrip- 
tion of a Fermi liquid based on the Landau theory. Since 
the Landau theory is valid only for small deviations of 
the quasiparticle distribution function n p from the Fermi 
function, only linearized version of the generalized hydro- 
dynamics can be constructed. Although such theory is 
formally equivalent to the linear response approach, the 
hydrodynamical formulation is generally more efficient 
for spatially inhomogeneous problems. In addition, the 
derivation of the generalized hydrodynamics of a Fermi 
liquid is of a general interest since it transparently shows 
how correlation effects contribute to the stress tensor. 

We_s±art from a kinetic equation for a charged Fermi 
liquidril 



d, 



'trip 



V£p dp 



where e p is the local quasiparticle energy, which depends 
on the distribution function and, consequently, is a func- 
tion of spatial coordinates. The gradient of this function 
gives an additional force which has its origin in the cor- 
relation effects. 

Following the derivation of the hydrodynamic equa- 
tions for a Fermi gas (Sec. II) we separate the macro- 
scopic and relative motion of a liquid by transformation 



if 



r(0) 



= TO* = TO 1 



1 



(58) 



dp 2 p ) V' 3" 

where F\ is the first of Landau Fermi-liquid parameters: 
m*pF f dfl' 



Fi 



47T 



fp> p Pi(cos6)i 



(59) 



where 9 is the angle between p and p', dfl' is an element 
of a solid angle around p' and Pi (cos 0) are Legendre 
polynomials. 

Zeroth and first moments of Eq. (|57]) correspond to 
the continuity equation and the equation for velocity of 
a Fermi liquid: 



d t Sn + n y i v l = 
mnodtVi + V jSPij — enoV.^ = 



(60) 
(61) 



In Eq. (61) we introduced the stress tensor of a Fermi 
liquid 



SPi: 



E 



PiPj 
TO* 



J TO 

P P,P 



\- , dn p 0) 
^ /p ' P d Pi 



vMp> (62) 



The first term in Eq. (|62j) corresponds to the kinetic 
stress tensor of quasiparticles, whereas the second term 
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describes the interaction of qusiparticles and is responsi- 
ble for a correlation contribution to the stress. 

According to the results of the previous section, the 
second order approximation (up to 7 2 ) gives a consis- 
tent hydrodynamical theory of a degenerate Fermi sys- 
tem. To make the system of equations (|60|)-(|6l|) complete 
one needs only an equation for the second moment or the 
stress tensor. To simplify the derivation of this equation 
we separate the stress tensor SPij into the scalar SP and 
the traceless %ij parts 



5Pa = SPSt-i + m- 



First we rewrite (p2|) as follows 



(63) 



(64) 



a quasiparticle. Thus the derivation of an equation for 
the stress tensor is straightforward. The final set of the 
generalized hydrodynamics equations of a Fermi liquid in 
the second order in parameter 7 ( ^6|) takes the form 



d t Sn + n V k v k = 
mnodtVi + ViSP + V ' jiTij 
d t 8P + KV k v k = 



enoVitp = 



(72) 
(73) 
(74) 

-v c iTij. (75) 



Hence, in a Fermi liquid we have again the system of 
equations which is similar to the Maxwell's theory of 
highly viscous fluids, but with the bulk K and the shear 
/i moduli directly related to microscopic Landau param- 
eters 



Ti 3 (p) =PiPj 



m pf 



77^ 



4?r Jp'pPiPj 



(65) 



In the Landau theory all momentums reside on the Fermi 
surface. By symmetry tensor Ty (p) is uniquely repre- 
sentable in the form 

T iS (p) = ap 2 5ij + Pijpipj - ip 2 % ) , (66) 

where the constants a and (3 are determined by the fol- 
lowing two equations 



—T 

P F 



3a = 1 



1 t jl 2 R 

-T-PiT ijP:j =a + ~(3- 

Pf 

1 



m pf 



7T Z 



m pf 



7T 



dZV_ 

4tt /p ' p ' 



(ppO 



'\2 



4tt /p ' p - ! 



Pf 



Solutions to these equations are 



a = 1 + F , 



= 1 + -F 2 
5 



(67) 



(68) 



(69) 



Substituting a and /3 ( |69| ) into Eqs. (p6|) , (|64|) we obtain 
the following microscopic expressions for SP and tt^ 



5P = (i + F )J2£-;K = v s 5 Pi 



(70) 
(71) 



wherej-uf = (1 + F^)p 2 F / imm* is the square of velocity of 
soundEj and Sp — raSn is a variation of the mass density. 
The equivalence in Eq. ( f70| ) shows that the scalar part 
of the stress tensor SP exactly equals to the variation of 
the pressure of the Fermi liquid. 

Expressions for SP and nij (|7(]) and ( |7l| ) have the same 
form as for the gas case except the factors in front of the 
sums and the change of the bare mass to the mass of 



2 5 1 
K = mn v s = g-Po — 



Fn 



1 J? ' 



1 + 

1 + 



(76) 
(77) 



where Pq is a pressure of an ideal Fermi gas. 

The bulk modulus ((76|) and the shear modulus (|77| ) 
exactly coincide with those of Ref. [l5| where they have 
been obtained by analyzing the linear response function. 
In ideal Fermi gas (Fi = 0) Eqs. ( |76| ) and ( |77j ) transform 
to the ideal gas moduli (see Sec. III). 



V. HYDRODYNAMIC THEORY OF 
COLLECTIVE MODES 

A. Plasma waves in a homogeneous system 

The frequencies of collective plasma modes in a 
charged liquid have an order of magnitude of oj p which 
commonly resides in the high frequency (collisionless) re- 
gion uj/v c 3> 1. The plasma oscillations are solutions 
of the linearized system of the hydrodynamics equations. 
To obtain the dispersion of plasmons we have to solve 
Eqs. ([72] )- ([75|) with the collision term being omitted in 
Eq. (|75|). I n the collisionless limit it is convenient to in- 
troduce a displacement vector u and rewrite this system 
in the form of the elasticity theory (]50|), ( |5l|) , where K 
and n are determined, in general, by Eqs. (76) and J7^). 
Substituting Eq. (51) to Eq.(^o|) we have an equation for 
the displacement vector of a "charged elastic medium" 



mnod t u — (K 



-M)V(Vu) 



en^Vip = 0, 

(78) 



which should be solved together with the Poisson equa- 
tion 



V ip = -47reVn u 



(79) 
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Considering a plane-wave solution e l ^ wt qr - ) we get the 
dispersion of the plasma waves 



where 



2 2 



4/i 



(80) 



(81) 



The first contribution in Vq equals to the square of the 
sound velocity (the result of the ordinary hydrody- 
namics) and comes from the fluctuations of pressure 5P 
(see formulas (|7^) and (|76|)). The second contribution 
arises from the traceless part of the stress tensor 7Ty. 
The sum gives a correct coefficient which can be obtained 
from the longitudinal response function of a charged 
Fermi liquidJlia In the case of a nonintcracting gas at 
arbitrary degeneracy the expression for Vq reduces to the 
formula 



3muo 3mno mno 



=<v; > 



which coincides with the result of the RPA (up to the 
second order of qu/u> ~ 7) and in a degenerate Fermi gas 
gives the well known result §ff-- It is straightforward to 
show that taking into account the third and fourth mo- 
ments in the linearized system of equations leads to the 
correct plasmon dispersion up to q A . 



B. Plasma oscillations of a confined electron liquid 

In the case of one dimensional motion (e.g. for plasma 
waves in an infinite medium which are considered in the 
previous subsection) the displacement vector u depends 
only on one coordinate. Hence, the equation of motion 
takes the following simple form: 



ySV u Wtp = 0. 

m 



(82) 



Eq. ( J82| ) coincides with that of the linearized Bloch's 
theory. The only difference is that the dispersion coeffi- 
cient equals to the correct value instead of a square 
of a velocity of sound v^. It seems easy to phenomeno- 
logically improve the BHT. One could try to replace the 
dispersion coefficient by the correct value. Qp-and hence 
obtain a correct dispersion of plasmon .ffl'LEil However, 
from the general point of view such an improvement can- 
not be consistent since the correct static limit requires 



vj. as a coefficient in Eq. (|S2|). The frequency-dependent 
coefficient was commonly used in hydrodynamics calcu- 
lations to recover both the high- and the low-frequency 
regimes (see ReL- fl] and references therein) . It has been 
then recognizedliiJthat the theory with a frequency de- 
pendent dispersion coefficient violates the harmonic po- 
tential theorem (HPT) in a spatially inhomogeneous case. 
Nonetheless, it was believed that since the theory with 



the replaced coefficient gives the correct plasmon disper- 
sion in the infinite medium it should hopefully give the 
proper description of the high-frequency plasma os^l^i- 



tions for any geometries of a confined electron gas. I 

As we show in the previous sections, the consistent 
theory which is valid both in a high- and low-frequency 
regimes must inevitably include the tensor fields of higher 
rank (e.g. 7Ty in the second order theory). At low fre- 
quency TTij goes to zero that gives a correct static limit. 
However, in the high-frequency regime it does not con- 
tain any small parameter and gives the contribution of 
the same order of magnitude as the variation of the scalar 
pressure. Within the linearized equations the contribu- 
tion of the tensor 7r.y can be interpreted as a shear stress 
of a Fermi system. Formally this leads to more compli- 
cated structure of the differential operator in Eq. (|78| ) 
than one of the operator expected from the BHT Eq. 

In the present subsection we consider two simple ex- 
amples of a collective motion of a Fermi liquid which is 
confined in three and one dimensions. These examples 
show that the collective eigenfrequencies strongly differ 
from predictions of the standard hydrodynamics. 

Consider a homogeneous (in equilibrium state) elec- 
tron liquid which is confined by an external potential. 
The homogeneity of the equilibrium state can be reached 
for example with the help of a positively charged back- 
ground of the density equal to the density of electrons 
uq or with the parabolic potential of a proper frequency 
loq. For example, in a three-dimensional spherically 
symmetric case (parabolic quantum dot) the frequency 
ujq = uip/ -\/3. The two-dimensional (parabolic quan- 
tum wire) and one-dimensional (quantum well) potentials 
have frequencies loq — uo p /^/2 and ujq — <jJ p respectively. 

Eigenmodes are solutions to equations j5fj|) , (|5l"l) , j79|) 
with the prop er boundary conditions. Integration of Eqs. 
( |50| ) and (f79|) across the surface gives the set of boundary 
conditions 



n\7ip>(r s ) - nV^ < (r s ) 



0, 

47ren nu(r s ), 



(83) 
(84) 
(85) 



where r s is a coordinate of a surface point, </? > and ip < 
are respectively the electrostatic potential ip outside and 
inside the surface and n is the unit vector normal to the 
surface. The boundary condition ( p3[ ) corresponds to the 
"free surface" and allows an electron liquid to cross the 
boundary and oscillate with respect to the equilibrium 
position. For example, the rigid oscillations of an elec- 
tron liquid u = Uoe~ iuJt (uq is a constant vector) satisfy 
Eqs. @),@,@) and conditions (||)-(||). Such oscil- 
lations have the frequency u> = loq. This fact apparently 
reflects the HPT. 

The above system of equations and boundary condi- 
tions also satisfy the HPT (or generalized Kohn theorem) 
in the most general formulation, which was suggested by 
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DobsonB We demonstrate this for a Fermi liquid con- 
fined by a spherically symmetric harmonic potential with 
luq = LOp/ \/3 and subjected to external homogeneous elec- 
tric E and magnetic H fields. The HPT motion corre- 
sponds to rigid displacement u(r, t) = uo(i) of the liquid. 
Vector uo(£) is obviously the center-of-mass coordinate. 
Since the stress tensor is equal to zero, the boundary 
condition (|8q) is trivially satisfied and equation of mo- 
tion (]50|) has the form 



md t u — eH x d t u + eE — eV( ) c < = 0. 



(86) 



The scalar potential ip satisfies the Laplace equation and 
electrostatic boundary conditions (|84|), fl85|). For rigid 
motion these conditions coincide with the conditions on 
the boundary of a homogeneously polarized sphere with 
a polarization vector P = — enoiio. The soliton of this 
problem inside the sphere is a depolarization field 



-Uq. 



r-, < 4 7T 
\7lf< = — P 



Substituting this solution into the Eq. (J8q) we find that 
the center-of-mass moves according to equation 

md 2 u - eH x <9 t u + eE + u 2 u = 0, (87) 

in exact correspondence with the HPT.0 

To avoid a confusion we would like to outline that the 
boundary condition ( |S3| ) assumes the parabolic potential 
outside the region occupied by an electron liquid. Physi- 
cally it can be realized, for example, at the boundary of a 
depletion region in p-n junctions or in artificial quantum 
structures. The condition ( j33| ) is not directly applica- 
ble to a model surface of metals with abrupt change of 
the positive background at the boundary. The reason is a 
jump of a potential at the surface which forbids electrons 
to penetrate the boundary. 

Let us separate longitudinal and transverse compo- 
nents of the displacement u 



All solutions of Eqs. 
ary conditions 

(cf - 2cf) iii V 2 ip + 2cjn i VjVi^ 
+ 2c?n i (V<t4 + V J -«j) 
nV(ip < + Anenoip) = nV(/? > , 



9fi-(B2h must fulfill the bound- 



(93) 
(94) 
(95) 



where all functions are taken at the surface. 
(p9|)-(|95|) we introduced the notations: 

<?i=v 2 = ——(K+-ii), 
mno 6 



mno 



In Eqs. 

(96) 
(97) 



The boundary condition (f53|) demonstrates a peculiar- 
ity of the generalized hydrodynamics. Longitudinal and 
transverse oscillations are, in general, mixed at the sur- 
face. This effect is well known in the elasticity theory!] If 
the shear modulus is nonzero, the incident purely longi- 
tudinal (transverse) wave transforms into the mixed wave 
containing both longitudinal and transverse components. 
As a result the structure of eigenmodes is changed in com- 
parison with the common hydrodynamics (/i = 0). The 
qualitative structure of the wave is not changed only in 
the case of a normal incidence (for example radial modes 
in the spherical symmetric system which are purely lon- 
gitudinal). However, even in this case the shear modulus 
contributes to the boundary condition. The boundary 
condition for a purely longitudinal motion takes the form 



(cf - 2c 2 )n t X/ 2 ip - 



2cfnjV j-V 'iip = 



which strongly differs from the common hydrodynamical 
condition V 2 ip = for a parabolic potcntialE 

Bellow we solve equations (89)-(^5|) for two specific 
cases and show that nonzero shear modulus strongly in- 
fluences the dispersion of the eigenmodes even for purely 
longitudinal waves. 



u = Vip + u t , Vu f = 



(88) 



("0 is a potential of the longitudinal c omp onent) and 
rewrite the system of equations ( |50"| ) , (|5l]) , ([79]) and bound- 
ary conditions (|S3"|)-(|S5"|) in terms of variables ip and u t . 
We are looking for harmonic in time solutions ~ e~ lut . 

In the region occupied by an electron liquid we have 
the system of differential equations: 



u 2 -0 + c?V 2 V> + — ¥> = 0, 



LU 2 U t 



c 2 V 2 u t = 0, 



Vu t = 0, 



V 2 {ip + 4xen ip) = 0. 



(89) 

(90) 
(91) 



Outside the surface the Laplace equation for the scalar 
potential ip must be satisfied 



V> = 0. 



(92) 



1. Radial plasma oscillations in a parabolic quantum dot. 

In a spherical symmetric quantum dot the normal vec- 
tor n = f . Thus all normal derivatives transform to the 
derivative with respect to the radial coordinate r. The 
direction of the displacement u for radial modes coin- 
cides with f . Hence the transverse component u t equals 
to zero and functions i\) and tp depend only on r. 

Radial solutions of Eqs. (f§), ©, @ which are reg- 
ular at zero and at the infinity take the following general 
form 



, , . . sm qr „ 
^{r)=A 1 + </?<(r) 



<P>(r) = -, 



A 2 S -^- + B 2 , (98a) 



(98b) 
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Substituting (98|) into Eqs. (p9|), ( |91|) and into boundary 
conditions for the scalar potential fl94|) , (|95|) one gets the 
result 



sin 



sin gi? 



r 



iy9< (r) = AiAireno 



uj 2 v + c 2 g 2 i? 

sin qr sin qR 
~r R~ 



(99) 
(100) 



Here R is the radius of the quantum dot. The wave vector 
q is related to the frequency u 



2 2,22 



(101) 



The last boundary condition (|93J) in the spherical sym- 
metric case reads 



2c?)V ! ty(i*) + 2c 2 d 2 .^{R) = 0. 



This equation together with the solution (|99|) gives the 
final dispersion equation 



tan qR 



qR 



1 - (q/2q) Vi? 2 



(102) 



which determine the allowed values of q and, conse- 
quently, the frequencies of the eigenmodes (101). The 
corresponding result of BHT can be obtained from Eq. 
( |102| ) in the limit ct — > and reads 

tan qR = 0. 

The quantity (ci/2ct) 2 which governs the difference of the 
last two dispersion equations can be expressed in terms 
of Landau parameters: 



/ /o \2 1 (, 5 1 + F 



This value is obviously far from infinity (as assumed 
in the BHT). For example, in the case of a Fermi gas 
(Fi = 0) we have (c;/2c t ) 2 = 3/4 which is less then one. 



2. Surface plasma modes at the edge of a parabolic potential 

well. 



Let us consider the surface plasma oscillations at the 
edge of an electron system confined by a one dimensional 
parabolic potential with loo = <jJ p . The situation is com- 
mon in p-n junctions where the potential in the deple- 
tion region is parabolic and density of electrons (holes) 
in n (p) regions is approximately constant up to the cor- 
responding boundary of the depletion region. Possible 
applications of the results to the surface modes of the 
parabolic potential well is obvious. Within BHT this 
problem has been considered in Ref. |^. 

We suppose that the electron gas occupies the lower 
half space (z < 0). Let the surface wave propagate 



along cc-axis ip,u t ~ e lqx . Due to the symmetry, the 
y-component of the vector u t equals to zero. Hence the 
solution of Eq. (Ml) can be taken as 



u t = ae KZ e lqx 



(103) 



where a = [a x , 0, a z ) is a two dimensional constant vector 
and 



2 2 2 /2 

K =q -u) /Ct 



(104) 



Due to the condition Vu t = (see Eq. (|90|)) the con- 
stants a x and a z are not independent: 



-><1 



The general solution of Es.(B9, WW) takes the form 



^ = {Ae qz + Be pz ) e tqx , 

| ^-^Ae qz + 4:nen Be pi 



(105) 
(106) 



with 



V 



q 2 + {ul-u 2 )/c 



In the upper half-space we have the solution of the Pois- 
son equation (pJ 



Lp> = Ce 



(107) 



The relationship of the constants B and C to the constant 
A can be easily obtained from the electrostatic boundary 
conditions (p4[), (p5|). The result for the potential of the 
irrotational part of displacement tp is 



ip = A 



1 



(108) 



where uj s = (J p /y2 is the classical surface plasmon fre- 
quency in the infinite wave length limit. Let us rewrite 
the boundary condition (^3|) explicitly in the case of the 
plane boundary : 

(c 2 -2c 2 KV 2 V> + 2c 2 d 2 V> + 4c 2 <9 z < = 0, 
c\(d x d z i\) + d x u z + d z u x t ) = 0. 



Substitution of the solutions (103) and (108) to the 
boundary conditions gives the system of equations for 
the frequencies of the surface plasma modes: 



1 — — 77 ) (cj 2 — lj 2 — 2c 2 g 2 ) A - 2c 2 t na z = 



2 2 



At q = we have 



(109) 

A + c 2 (k 2 + q 2 ) a z = 0. (110) 



uj 2 

u 2 =u 2 s = y> az = 0, 
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which corresponds to the usual infinite wave length sur- 
face plasma oscillations. At small q the dispersion of the 
surface plasmon takes the form 



U3 



2c 2 q 2 



(112) 



0J S 



2c 2 t q 2 + i 



Ac 3 q 3 



(111) 



The damping of the surface plasmon is not surprising 
since at q ^ there is mixing with the transverse 
modes which can not decay from the surface at high 
frequency. There are decaying transverse solutions only 
for ui 2 < c 2 q 2 . Hence only the propagating transverse 
waves are allowed at the frequency ~ w a . The coupling of 
the longitudinal and transverse components leads to the 
energy transfer to propagating waves and, consequently, 
provides damping of the surface plasmon. It should be 
mentioned that the hydrodynamic theory does not take 
into account the Landau damping due to the single par- 
ticle excitations. Inclusion of the Landau damping will 
bring the usual contribution to the lifetime of the sur- 
face plasmon. Yet in a weakly coupled electron gas the 
Landau damping destroys the transverse modes. More 
precisely, if Ct < vf the transverse modes are damped. 
Hence the decay channel of the surface plasmon into the 
transverse waves will be really open if the inverse inequal- 
ity ct > vf is fulfilled and transverse modes are well de- 
fined. Using Eq. (^) for Cj and Eq. ( f77| ) for /i we can 
express the condition ct > vp in terms of Landau Fermi 
liquid parameters: 

F 1 (l + F 2 /5) + 3F 2 /5 > 6. 

Consequently, the transverse modes and the decay chan- 
nel described above could exist for strongly correlated 
charged Fermi liquid. It was mentioned in Ref. |l5| that 
the condition should be satisfied for a electron liquid of 
low density. 

With the accuracy up to q 2 w e ca n neglect the imagi- 
nary term in the dispersion law (111). In this long wave- 
length approximation the surface plasmon is a longitudi- 
nal wave with the dispersion 

I 

D t P + \PV ■ v + hit ■ V) • v + Tr(V • L (3) ) = 0, 



The last equation shows that the dispersion of the al- 
most longitudinal surface plasmon is totally determined 
by the shear modulus (c 2 ~ fi). Within the BHT the 
shear modulus equals to zero. In the case c* = equa- 
tions (109, 110) lead to the dispersionless surface plasmon 
with the frequency uj — lu s . This agrees with the results 
obtained inEI for a wi de p arabolic quantum well. The cor- 
rect dispersion law (112) obtained from the generalized 
hydrodynamics has the dispersion coefficient 2c 2 which is 
of the same order of magnitude as the coefficient v 2 = c\ 



in the volume plasmon dispersion (|30|). For example, in 
the case of a high density electron gas (Fi » 0) we have 
the ratio 



2c 2 /c 2 



VI. 



THE FOURTH-ORDER THEORY. THERMAL 
CONDUCTIVITY 



As we have seen in Sec. Ill, the theory of the second 
order in parameter 7 ( p6| ) does not contain effects of ther- 
mal conductivity which are responsible for the thermal 
equilibration processes. Hence the second-order approx- 
imation cannot give correct static solution which should 
correspond to a constant temperature (VT(r) = 0). To 
account for the heat conduction one has to consider at 
least fourth-order approximation which contains two ad- 
ditional equations for the third L^ 3 ) and fourth L^ 4 ) mo- 
ments. The system of equations follows Eqs. (p2])-(p5|). 
The first two equations - the continuity equation and the 
equation for the velocity field, are the same as in the 
second-order theory (Eqs. (j27|) and (|28| ) respectively). 
In the equations for the scalar P and traceless 7r parts 
of the second moment we get additional terms which are 
proportional to the gradient of the L^ 3 ^ (see Eq. EJ): 



(113) 



D t -K + tt(V • v) + {(tt • V) <g> v} s - l|(7T ■ V) • v + P ({ V <g> v} s - l^V ■ v ) + V ■ L (3) - Tr(V ■ L (3) ) = I (2) . (114) 



To close the set of equations with the accuracy of 7 4 we should take an equation for L^ 3 ) in the general form (|25|) and 
omit the spatial derivatives of the fifth and the second moments in the equation for L^ 4 ': 



AL (3) + L< 3 )(V • v) + {(L< 3 ) ■ V) ® v} s + V • L< 4 > 
AL (4) + L< 4 > (V • v) + { (L^ ■ V) ® v} s - iW . 



— {l( 2 )®(v-l( 2 >)} s =i( 3 \ 



(115) 
(116) 
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Eqs. (g7|), and (|113|)-(|116|) constitute the system of 
the generalized hydrodynamics equations in the fourth- 
order approximation. 

The last term in the left-hand side of Eq. flll3| ) can be 
rewritten as a divergency of a vector 



Tr(y-L^)=V J L^J=V J Q J , 



(117) 



(3) 

where the vector Qj = describes the flow of energy 
in a comoving frame. 

To demonstrate the recovering of the thermal conduc- 
tivity we take the moments of the collision integral I^ 3 ^ 
and I^ 4 -* in a linear approximation, similar to (|3q): 



l( 3 )=-, cl L( 3 \ 



(118) 
(119) 



where we introduced the notation hp for the fourth mo- 
ment of the distribution function f£(r,t) (see Eqs. (|3l|)- 
(ID): 



(120) 



Due to the spherical symmetry of the function f£ the 



fourth rank tensor Lj? (120) takes the form: 
Lf = {SijSki + SkjSn + SijSik)LF, 



(121) 



Parameters v c \ and v C 2 in Eqs. (|118| ) and ( 119 ) can be 
considered as phenomenological coefficients which are, in 
general, not equal to the parameter v c d36|), but have the 
same order of magnitude. 

In t he l ow frequency limit lo/v c « 1 we solve Eqs. 
(114)-( |TT6| ) using a perturbation theory. Since the colli- 
sion terms in all these equations have an order of magni- 
tude of the collision frequency the solution in zero order 
of l/v c follows from the equations 



lO)=0, j= 2,3,4 
and takes the local equilibrium form 

= L F , =0, 7T = 0. 



The stress tensor becomes diagonal and proportional to 
the local pressure 

I/ 2 ' = 1P F 



The first order solution of equation (30) for the tensor 
7r has been already obtained in Sec. III. It has the form 
of the viscosity tensor (f4l|). Eq. (114), which determines 
7r with the accuracy of 7^, differs from Eq. (^) since the 
thi rd ra nk tensor L/ 3 ) contributes in the left-hand side 
of (114). H owev er, the contribution of the third moment 
in the Eq. (114), being of the higher order in oj /v c , does 
not change the first order result of Sec. III. Hence the 
viscous solution ([y]) is still valid. 

Let us consider solution of Eqs. (115), ( |116| ) for the 
third moment L/ 3 ). In the first-order of u. 1 jv c only the 
last two terms in the left-hand side in Eq. (115) con- 
tribute to the equation 

- v cl \}^ « V ■ l/ 4 > - — {l/ 2 > ® (V • L( 2 ))}s, (122) 
mn 

where zero-order expressions for L/ 4 ) and L/ 2 ) should be 
used in the right-hand side. 

Substituting zero-order (local equilibri um) solutions 
(L< 2 ' = 1P F and L< 4 ' = L F ) to Eq. |L2§) one gets 
the the first-order expression for the third moment 



.(3) 



-— 1 1 <g> (vL F - —P F VP F 
v c \ I V mn 



(123) 



Contraction of the tensor L/ 3 - 1 (|l23j) over the couple of 



(3) 

indexes gives the energy flow vector Qj = L\-l, which 



enters the equation of the energy conservation ( 113 ) 

(124) 



5 / I 
Q = VL F P f VPf 

v c i \ mn 



Let us show that vector Q is proportional to the gra- 
dient of temperature VT(r, t) at arbitrary degeneracy of 
an electron gas. The functions Lf and Pf depend on 
spatial coordinates via the spatial dependence of the lo- 
cal chemical potential n(r 7 1) and the local temperature 
T(r,t) entering the distribution function /p(r,i) 



ffct) = f 1 



T 



P 

2m 



Hence gradients of Pp and Lp can be rewritten as follows 



3m 
1 



VT 1 



T 6m 2 



2^ p -*r 



T J 15m 3 



de r 



VT 1 

T 30m 4 



(125) 
(126) 
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The partial integration a llows one to express the mo- 
mentum integrals in Eqs. ( |l25| ), (|l26| ) in terms of macro- 
scopic variables n, Pp and Lp: 



3m ^— ' *" cfe 
p 



■E/. 



E^ = -5-E^/p = -15- 2 ^, 



30m 



—y 



0s 
A 1 



9e D 30m 3 
p F p 



7 E^ = -K 



Using the last expressions we can represent Eqs. (125) 
and (126) in the following compact form 



VP F 



V 1 I m 2i 



(127) 
(128) 



Comparison of Eqs. (|127|) and (|12§) with Eq. QlgJ ) 
shows that the contribution of the gradient of the chem- 
ical potential exactly cancels in the equation for the en- 
ergy flow vector. Thus the vector of the energy flow takes 
the usual form 



Q = -AVP 
with the coefficient of heat conduction 



A 



2^ c iT 



7L F - 



5P 



(129) 



(130) 



Finally Eq. (pT3| ) with tt @ and Q ( po| ) transforms 
to the common hydrodynamical equation of the energy 
conservations 

DtP + ^PV-v + H(7r.V)-v + V-Q = 0, (131) 

which contain both the viscous and the heat conduction 
terms. 

The formula fll3C| ) for the coefficient of heat conduc- 
tion is valid at arbitrary degree of degeneracy of an 
electron gas. In the nondegenerate case Pf = nT and 
Lp = nT 2 /m. Hence Eq. (130) transforms to the result: 



A = 



5nT 



In a degenerate gas at T = we have Lp — hPp/lmn 
hence the expression in the brackets in ( 130 ) cancels. At 
T/sp <C 1 the correction has the order of ma gnitu de 
(T/ep) 2 and the coefficient of heat conduction ( 130 ) is 
proportional to T jv c \. Since the collision frequency in a 
degenerate Fermi system is a quadratic function of T the 
heat conduction diverges as l/T (see, for example, Ref. 
1). 



Thus, in the low frequency li mit t h e gen eralized hydro- 
dynamics equations (|27|), ( Pq , (113)- (116) coincides with 
the correct set of the common hydrodynamics which con- 
sist of the continuity equation, Navier-Stokes equation 
and the equation of the energy conservation in the form 
(11). 

In the high frequency (collisionless) limit the general- 
ized hydrodynamics leads to the following plasmon dis- 
persion 



< v ;> 



< V l > 2 4 
n 1 



V 



that is exactly the result of the linear re spon s e the ory. 

The system of equations <M) , (|rn|)-(|TT6|) with 

the collision terms (|6|), (118), (|119[) constitute the closed 



set of equations and provides the smooth interpolation 
between high and low frequency regimes. 



VII. CONCLUSION 

A description of dynamics of many-electron systems in 
terms of macroscopic collective variables, which is usu- 
ally referred to as a hydrodynamical approach, provides 
a simple, physically transparent and powerful tool for 
studying spatially inhomogeneous problems. The com- 
mon Bloch's hydrodynamics of an electron gas is based on 
the assumption of a local equilibrium and consequently 
shows a number of inconsistencies. We have shown that 
it is possible to construct an inherently consistent gener- 
alized hydrodynamics which correctly describes both the 
collisionless high-frequency limit and the collision domi- 
nated low- frequency regime, where the theory coincides 
with the standard Navier-Stokes hydrodynamics. The 
theory follows from the long-wavelength expansion of a 
kinetic equation and requires inclusion of new collective 
variables with a nontrivial tensor structure. We remind 
that only scalar (pressure and density) and vector (ve- 
locity) variables enter the common hydrodynamics. 

The appearance of higher rank tensors is physically 
evident since the occupied region in a momentum space 
loses its spherical symmetry under a general (not locally 
equilibrium) evolution of the system. A need to describe 
a non-spherical isoenergetic surface inevitably requires 
tensor variables. As long as this surface is smooth it is 
possible to approximately describe it as a surface of a fi- 
nite order and thus to take into account the tensor fields 
of the finite rank. We have actually shown above that 
the smoothness of the surface, which bounds the occu- 
pied region in the momentum space, is governed by the 
basic parameter (26). 

In extension of the previous publication,!!!] we pre- 
sented a new general and more transparent derivation 
of the generalized hydrodynamics and showed how the 
closed set of equations of the standard Navier-Stokes hy- 
drodynamics (including the heat transport) is recovered 
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in the low frequency regime. We have shown that the gen- 
eralized hydrodynamics can be built on the basis of Lan- 
dau theory of Fermi liquid which allows to determine the 
correlation contribution to the stress tensor. Although 
this hydrodynamics (as the Landau theory of Fermi liq- 
uid itself) is meaningful only in the linear approximation, 
we belive that it could be of a practical importance since 
it allows to express the results of hydrodynamical cal- 
culations, which are relatively simple, in terms of the 
microscopic Landau parameters. For example, the the- 
ory offers a possibility to determine the contribution of 
Fermi-liquid correlation effects to the eigenmodes of spa- 
tially inhomogeneous systems. 

In a linear approximation the tensor variables, which 
describe the absence of a local equilibrium, can be inter- 
preted as an effective shear stress of a liquid. On a time 
scale much longer than the collision time the contribu- 
tion of the shear stress vanishes and the dynamics of a 
liquid is governed only by the usual bulk stress. How- 
ever, as the collective modes of a charged liquid normally 
belong to the high frequency range, the shear contribu- 
tion cannot be neglected. It is important to realize that 
this contribution cannot be modeled, even qualitatively, 
merely by a change of the bulk modulus - the procedure, 
which has been frequently used in, literature to obtain 
the correct dispersion of plasmomaBLl The reason is that 
the relative contribution of the shear and bulk stresses is 
different for different modes and different geometries. As 
a result the eigenfrequencies and the structure of plasma 
modes in confined systems strongly deviate from predic- 
tions of BHT even with "improved" bulk modulus, which 
is adjusted to provide a correct plasmon dispersion in a 
homogeneous situation. To demonstrate this qualitative 
non-local-equilibrium effect we calculated eigenmodes 
of a Fermi liquid confined by a harmonic potential of 
different dimensionality. The results show a nontrivial 
contribution of the shear modulus to eigenfrequencies 
of the longitudinal plasma oscillation. The most trans- 
parent result concerns surface waves at the edge of ID 
harmonically trapped system. In the absence of the 
shear modulus these modes are absolutely dispersionless 
in agreement with the BHT results.Q They acquire a 
dispersion of the same order of magnitude as the bulk 
plasmon due to deviation from the local equilibrium. We 
have also found that existence of a nonzero shear mod- 
ulus may open a new channel of decay of these surface 
modes, in addition to the usual Landau damping. 
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